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ABSTRACT 


Nonlinear tapered flexbeam laminates, with significant ply waviness, were cut from a full- 
size composite rotor hub flexbeam. The specimens were tested under combined axial 
tension and cyclic bending loads. All of the specimens had wavy plies through the center 
and near the surfaces (termed marcelled areas), although for some of the specimens the 
surface marcels were very obvious, and for others they were much smaller. The 
specimens failed by first developing cracks through the marcels at the surfaces, and then 
delaminations grew from those cracks, in both directions. Delamination failure occurred in 
these specimens at significantly shorter fatigue lives than similar specimens without 
waviness, tested in ref. 2. 

A 2D finite element model was developed which closely approximated the flexbeam 
geometry, boundary conditions, and loading. In addition, the FE model duplicated the 
waviness observed in one of the test specimens. The model was analyzed using a 
geometrically nonlinear FE code. Modifications were made to the original model to reduce 
the amplitude of the marcels near the surfaces. The analysis was repeated for each 
modification. Comparisons of the interlaminar normal stresses, a n , in the various models 
showed that under combined axial-tension and cyclic-bending loading, for marcels of the 
same aspect ratio, a n stresses increased as the distance along the taper, from thick to thin 
end, increased. For marcels of the same aspect ratio and at the same x-location along the 
taper, a n stresses decreased as the distance from the surface into the flexbeam interior 
increased. A technique was presented for determining the smallest acceptable marcel 
aspect ratio at various locations in the flexbeam. 

NOMENCLATURE 


Eii, ^22 Young’s moduli in the 1- and 2- directions, GPa 
G 12 shear modulus, GPa 

L period of one marcel, mm 

N number of loading cycles 

P axial tension load, kN 

t(x) flexbeam half-thickness at distance x from fixed end, mm 

V transverse bending load, kN 

v transverse displacement at tip of flexbeam, mm 

y vertical distance from flexbeam surface, mm 

a waviness half-amplitude, mm 

6 transverse stroke of ATB test machine, mm 

8 flexbeam surface strain 

£ max maximum cyclic surface strain 

o n interlaminar normal stress in FE model, MPa 

Vi 2 Poisson’s ratio 


INTRODUCTION 


Tapered laminated composite flexbeams are used in helicopter rotor hubs to reduce 
weight, drag, and the number of parts in the hub. The taper in these flexbeams is achieved 
by terminating internal plies along the length. However, the manufacturing procedures 



required by these flexbeams can sometimes cause significant ply waviness throughout the 
tapered section of the flexbeam. 

In this paper the term marcel is used to describe an area where there is a sudden 
change in the ply direction, similar to a wrinkle in the ply. Because the manufacturing 
processes cannot always completely eliminate these marcels, the effect of this waviness on 
the flexbeam durability needs to be determined. Also, such a study on marcelling in 
flexbeams may be used to develop accept/reject criteria for these flexbeams. 

In ref. 1 , the effect of ply waviness on pull-off loads in composite hat stringer specimens 
was studied using the finite element method (FEM). Waviness was simulated in the FE 
model by an incline in the ply at the initial delamination location. Five different ply waviness 
angles were studied and the results showed that even mild amounts of waviness could 
result in significant reduction in stringer pull-off loads. 

Reference 2 described a method for determining the fatigue life of tapered composite 
flexbeam laminates. In that study, 25.4 mm (1 inch) wide coupon specimens were cut from 
a full-size flexbeam of S2/E7T1 glass/epoxy. The specimens were tested under 
combined axial tension and transverse bending loading in a servo-hydraulic load frame, 
called the Axial-Tension Bending (ATB) machine [2], using a frequency of 3 Hz and fully- 
reversed loading (R=-1). Under fatigue loading, delaminations typically initiated in the areas 
around the ply-drop locations and grew in both directions along the length. Typical fatigue 
lives were 10 5 to 10 7 cycles. 

In the current study, flexbeams of the same geometry, layup, and material, but with 
significant ply-waviness, were tested in an identical manner to those in ref. 2. The location of 
the delamination initiation and the fatigue lives of the wavy flexbeams were studied. A 
parametric study was also conducted using a finite element (FE) model of the tapered 
laminate to examine the effect of the severity of the ply waviness at a given location, and 
the effect of marcels of the same severity at different locations along the taper or through- 
the-thickness. 


EXPERIMENTS 


Specimen Configuration 

Tests specimens for this study were cut from a full-scale test flexbeam, which was 
manufactured in a closed-cavity tool. The full-scale test flexbeam and the cut coupon 
specimens are shown in fig. 1 . The full-scale flexbeam was cut in half cross-wise and then 
each section was cut lengthwise into 25.4 mm-wide coupons, yielding 6 test specimens. 
The coupons were symmetric with respect to layup and geometry and were designed with 
a nonlinear taper. The total number of plies in the laminates varied from 145 at the thick end, 
to 41 plies at the thin end. An idealized view (no waviness) of the cross-section of the 
upper half of the flexbeam is shown in fig. 2. This figure shows that the laminate had a 
woven fabric layer on the surface, four continuous “belt” sections, each 4 plies thick, and 4 
“dropped-ply” groups, on each side of the midplane. At the midplane was a symmetric 
ply-group, which was 5 plies thick. Each dropped-ply group had a maximum thickness of 
13 plies at the thick-end. The dropped plies were arranged in a non-uniform, staggered 
manner. The woven fabric on the surface was E-glass/E7T1-2 and the interior plies were of 
S2/E7T1 tape. Properties for both materials are given in Table 1 . The layups and material 
properties for each ply group shown in the figure are given in Table 2. 

Although all six specimens were cut from the same full-size flexbeam, the degree of 
waviness, and locations of the wavy areas, varied from specimen to specimen. 

Specimens 1-3, cut from one side of the flexbeam, had considerably more waviness than 
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specimens 4-6, cut from the other side. The ply waviness appeared as isolated areas of 
marcelling, and as waviness in the ply-groups at the center of the flexbeam. The waviness 
in the center ply-groups extended through most of the tapered region of the laminate, with 
the amplitude of the waves highest near the thick end, and gradually decreasing along the 
length. None of the specimens showed any waviness in the area of the laminate beyond 
the tip of Dropped-group 2 in fig. 2. An example of both types of wavy areas is shown in 
an edge view of specimen 1 , in fig. 3. The figure shows the waviness in the center ply- 
groups, and the enlarged photo shows two isolated areas of large amplitude marcels, one 
near each surface. The amount and type of waviness were fairly consistent on both sides 
of a given specimen, except for the large marcelled areas near the surfaces, which did not 
always appear on both sides. 

Specimens 1 -3, considered “severe” waviness, had waviness through the center-ply 
groups, as well as isolated areas near the surfaces where very large (high amplitude) 
marcels existed. Specimens 4-6, considered “moderate” waviness, had wavy plies 
through the center, but the waviness was of a lower amplitude than for specimens 1-3. 

Also, specimens 4-6 did not show the isolated areas of large marcels near the surfaces. 

Figure 4 shows an edge view of specimen 2, which had several isolated areas of 
extreme waviness near the surface, shown in the enlarged photos, in addition to the 
waviness through the center ply-groups. (The circular white marks near the surfaces in the 
photos are paint marks that were used to position the strain gages.) Since specimen 2 had 
the most severe waviness of the six specimens, and since the waviness, including the 
isolated surface marcels, was identical on both sides, specimen 2 was chosen as the basis 
for the finite element model. Specimen 3 was very similar to specimen 2, but did not have 
the marcel on the bottom surface in fig. 4. 

Figure 5 shows the edge of specimen 4. As the figure shows, the only evidence of ply 
waviness in this specimen was in the center ply-groups near the thick end. Specimens 5 
and 6 were similar to specimen 4. 

Axial Tension Bending Machine 

Specimens were tested under combined axial-tension and transverse-bending loading 
in the Axial-Tension Bending (ATB) Machine. The ATB, shown in fig. 6, is a servo- 
hydraulic load frame, which produces combined tension-bending loading. As fig. 6 shows, 
the axial load cell is located above the top grip, but below the pivot connecting the axial 
and transverse actuators. This allows the tension load to rotate with the specimen as the 
transverse load is applied. Flence, under axial load control, the magnitude of the tension 
load, P, remains constant as the specimen rotates under the transverse-bending 
displacement, 6. By controlling the axial tension actuator under load control and the 
transverse bending actuator under stroke control, a constant membrane load should be 
maintained throughout the loading cycle. 

Static Tests 

Initial static testing was conducted to determine the relationship between applied loads 
and the specimen deflection and surface strains. Specimens were first instrumented with 
strain gages at five locations along the length on each side: one near the junction of the thick 
and tapered regions (1 and 6), three along the tapered region (2-4 and 7-9), and one in the 
thin section (5 and 10). Figure 7 shows a schematic of the ATB and flexbeam, rotated 90° 
clockwise. The figure shows the numbered gages on each surface of the flexbeam and 


3 



Table 3 lists their location in distance from the fixed end. The specimen was clamped in the 
grips with the thick end in the fixed bottom grip. The gage length between the grips was 
165 mm (6.5 inches), and the specimen was placed in the grips so that within the gage 
section there was approximately a 12.7 mm (0.5 inch) thick region, a 127 mm (5 inch) 
tapered region and a 25.4 mm (1 inch) thin region, as shown in fig. 7. Note also in fig. 7, that 
the transverse bending load, V, was not applied at the top of the flexbeam, but was 
applied at the pivot point, which was 172 mm above the top grip. 

For static excursion tests, a constant axial tension load, P, of approximately 35.6 kN 
(8000 lbs.) was applied first. Then the bending load, V, was applied, in steps, to produce 
a transverse stroke, 6, in increments of approximately 2.54 mm (0.1 inch), up to a 
maximum stroke of 25.4 mm (1 .0 inch). At each transverse load step, the surface strains 
were recorded, as well as the transverse flexbeam tip-displacement, v. In order to 
measure the flexbeam tip-displacement, a spring-loaded direct-current differential tranformer 
(DCDT) was mounted to the side of the load frame. The DCDT detected the 
displacement of a bracket attached to the centerline of the top grip. 

In ref. 2, the peak surface strains in an identical flexbeam without waviness were plotted 
as a function of applied transverse stroke, 6. The relationship was shown to be linear as 
long as the axial load is held constant. Similarly, for this study, the peak surface strains were 
measured at gage 4 in fig. 7 (X=98.5 mm). The maximum surface strains are plotted in fig. 
8, as a function of the transverse load, with a constant axial load of P=8000 lbs. The results 
are shown to be linear for both specimens 2 and 4 (with severe and moderate waviness, 
respectively). 

Fatigue Tests 

Since the boundary conditions of the ATB differ from those of the full-scale flexbeam in 
the hub, it is more logical to control the fatigue tests to a desired maximum surface strain 
level, rather than a prescribed transverse deflection. Flence, for each specimen, a 
maximum strain level for fatigue testing was chosen, and then results of the type shown in 
fig. 8 were used to select the maximum cyclic transverse stroke, 6 (see fig. 7), to apply 
corresponding to the chosen strain level. The specimens from ref. 2, without wavy plies, 
were tested at transverse displacements of 27.9 mm or 30.5 mm, corresponding to 
maximum surface strains of 0.01 and 0.015 microstrain, respectively. However, because 
the static testing showed that the wavy specimens were susceptible to delamination at 
lower strain levels, the maximum strain levels (or maximum cyclic stroke) chosen for cyclic 
testing of these specimens were somewhat lower than those used for the specimens 
without ply waviness in ref. 2. Specimens 2-6 were fatigue tested at maximum strain 
levels of 0.0075 to 0.01 microstrain. 

In order to make the delamination damage easier to see, specimen edges were 
coated prior to testing with a thin layer of white paint. A constant tension load of 35.6 kN 
(8000 lbs.) was applied to the specimen, corresponding to the net axial stress due to the 
centrifugal force experienced by the full-scale flexbeam. The maximum cyclic transverse 
load, V, was applied by cycling sinusoidally to the desired maximum transverse stroke, at a 
frequency of 3 Hz, and using fully-reversed loading (R=-1). The specimens were cycled 
until they had extensive delamination damage along the length of the tapered region at one 
or more locations. Test results are discussed later in this report. 
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ANALYSIS 


Finite Element Model 

In order to duplicate the exact geometry of the wavy flexbeams, a software package 
known as MEGS (Modeling Exact Geometry from Scanned images) was used [4], With 
this software, a scanned image of the object to be modeled was used to create a wireframe 
image, incorporating as many “as manufactured” details as is desired. The wireframe data 
was then imported into a PATRAN file, to be used as the basis for developing a finite 
element model of this configuration. Because specimen 2 appeared to have the most 
severe waviness at locations near the flexbeam surface, as well as wavy plies in the center 
of the laminate, specimen 2 was chosen for use in the analysis and FE model. 

It was also necessary to apply the loads and boundary conditions to the model in a 
manner that duplicates the configuration of the ATB. A schematic of the configuration to be 
modeled by FE is shown in fig. 9. Fixed conditions were applied at the thick end of the 
composite flexbeam. Beyond the thin end of the flexbeam, additional elements were 
created to represent the upper grip and steel fixture connecting to the pivot point where the 
transverse load is applied (see fig. 7). The tension and bending loads are applied at the 
end of the model, rather than at the flexbeam tip. In this way, the loading conditions of the 
model duplicate the test conditions. Because the bending stiffness of the steel fixtures was 
two orders of magnitude greater than the flexbeam, the elements at the thin-end that 
represent the steel loading fixtures were modeled with a rectangular cross-section equal to 
the thin end of the composite flexbeam. A modulus was chosen such that the bending 

stiffness, El, was equivalent to the bending stiffness of the actual ATB fixtures of 1 .74x1 0 6 
N-m 2 (6.05x1 0 8 lb-in 2 ). 

Finite Element Mesh and Boundary Conditions 

A 2D finite element model of specimen 2 is shown in fig. 10. The model had a total of 
21 ,903 nodes and 6971 elements in the mesh. Eight-noded quadrilateral and 6-noded 
triangular plane-strain elements were used. A fine mesh, using one element per ply- 
thickness, was used for the ply-groups closest to the surfaces, in the areas around the 
marcels nearest to the top and bottom surfaces. A coarser mesh was used in the interior of 
the model and in the thinner region of the flexbeam. In those areas, ply-groups were 
modeled, rather than individual plies, and smeared properties were used. The smeared 
moduli in the global X-Y coordinate system are presented in Table 2. 

In order to assign appropriate material properties to the wavy plies, a local coordinate 
system was defined for each element in the model, with the 1 -direction parallel to the 
element side from the local node i to local node i+1 , as shown in fig. 10. The local t-n 
coordinate system was then used to define the material properties of each element. Figure 
1 1 shows enlargements of three marcelled areas near the surface. These have been 
designated marcels 1-3, as shown. The shaded areas in the enlarged figures represent 
resin-rich regions. The elements in the resin-rich areas (not shown in fig. 11) were assigned 
neat resin properties, as listed in Table 1 . 

The u- and v-displacements at the nodes at the thick end of the model were prescribed 
zero values to simulate clamped-end conditions. An axial tension load of 35.6 kN (8000 
lbs.) was applied at the thin end (X=337 mm) as a concentrated load (See fig. 7). 
Transverse bending was produced by a point load of V=-4.45 kN 
(-1000 lbs.) applied at the thin end of the model, corresponding to the pivot point in the 
ATB load frame. To determine the effect of the bending load on each surface, the bending 
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load was applied first in the negative Y-direction, and then in the positive 
Y-direction. 

Computational Methods 

The ABAQUS finite element code was used in the analysis. Because the flexbeam 
undergoes large deflections, the geometric nonlinear solution option was used. Also, as 
with the ATB load frame, the axial load in the model was able to rotate with the flexbeam as 
it deformed under the transverse load. The ABAQUS program was used to calculate 
displacements, and internal stresses and strains. 

To study the effects of changes in the ply waviness, modifications were made to the 
original mesh. Using the PATRAN code, marcels 1 , 2, and 3 were modified several times 
to reduce the amount of waviness at each location and the resin-rich areas were removed. 
The new models were analyzed with the same loading and boundary conditions, and the 
results were compared. 


RESULTS AND DISCUSSION 
Global Response Comparison 


Calculated values from the FE model for transverse displacements, v, at the 
flexbeam tip were compared to the test results to determine the accuracy of the FE model 
to reproduced the global behavior of the test specimens under the same loading 
conditions. Figure 12 compares the maximum surface strains vs. flexbeam tip 
displacement, as measured by gages 4 and 9, in the static excursion tests, along with the 
calculated results from the ABAQUS FE model. The transverse displacement was varied 
while the axial load was held constant at 35.6 kN (8000 lbs). The agreement is good 
throughout the range of tip-displacement, on both the tension and compression sides, 
although the calculated strains are slightly lower everywhere. 

In fig. 13, measured surface strains from the ten strain-gage locations are compared to 
the ABAQUS calculated strains along the flexbeam length. The calculated values are 
shown in the solid circles and test results are shown for a severely marcelled specimen 
(specimen 2, open circles) and a moderately marcelled specimen (specimen 4, open 
squares). The results shown are at P=35.6 kN (8000 lbs.) and 6=27.9 mm (1.1 inch). 
Gages 1-5 were on the tension surface, and gages 6-10 were on the compression surface, 
as shown in the inset in fig. 13. The agreement is reasonable for gages 1 , 2, 4, and 5 on 
the tension surface. Flowever, at the gage 3 location (X=85.8 mm, 3.38 in.), the calculated 
strain is much higher than the measured strain from specimen 4. The FE model duplicated 
specimen 2, which has a large marcelled area near the surface at gage 3, which is not 
present in specimen 4. Flowever, the data point for specimen 2 at the gage 3 location is 
shown with an upward arrow to indicate that the strain reading had exceeded the limit of the 
gage (.012 mm/mm). As the figure shows, the surface strains varied for the different 
specimens. Along the compression surface, FE results are not as accurate near X=0, but 
are reasonably close at gage locations 8-1 0 and follow the trend of the data. Based on 
these results, the FE model appears to duplicate well the global response of the test 
specimens under loading in the ATB. 
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Static Test Results 

During the static excursion tests of specimens 1 , 2, and 3, all of which had severe 
waviness, some delamination damage occurred. This was first noticed as faint “clicking” 
noises as the transverse bending loads were applied. A visual check was made before 
increasing the load each time. In specimen 1 , a crack formed first at the location 
corresponding to marcel 1 in fig. 1 1 , at a transverse displacement of 6=12.7 mm (0.5 in). 
Several small delaminations were observed through-the-thickness at the same location. As 
static loading continued, more delaminations developed at other interfaces through-the- 
thickness. Figure 14 shows a photo of the coated surface with the delamination damage. 
Along with the internal delamination damage, there was significant surface ply splitting in the 
surface fabric above the surface marcel. Because of the large amount of delamination 
sustained in the static testing, specimen 1 was not fatigue tested. Specimens 2 and 3 each 
developed a small delamination at marcel 1 also, approximately one quarter of the 
thickness from the surface, when the transverse displacement reached 6=20.3 mm (0.8 in.). 
These specimens were fatigue tested with the initial delamination damage. Specimens 4-6 
had no apparent damage after static testing to a maximum transverse displacement of 
6=27.9 mm (1.1 inch). 

Fatigue Test Results 

The specimens were visually monitored throughout the fatigue loading cycle. For 
specimens 2 and 3, which had a large marcel at the location “marcel 1” (see fig. 11), 
damage began at that location, as a crack through the wavy area, followed by a 
delamination at the interface under the marcel. This initial crack formed almost immediately, i. 
e., at less than N=50 cycles. As cycling continued, delaminations formed at neighboring 
interfaces, through the wavy areas. These delaminations grew in both directions with further 
loading. The internal delamination damage was always accompanied by splitting and 
peeling of the surface fabric ply at the location over the surface marcel. Continued fatigue 
loading caused delaminations to form and grow at interfaces in the wavy center-ply-groups, 
near the thick end of the flexbeam. Figure 15 shows the delamination damage in 
specimens 2 and 3. 

The first observed damage in specimens 4-6, which initially did not appear to have any 
marcels near the surfaces, appeared as a crack near the surface, similar to the initial damage 
in specimens 2 and 3. These cracks occurred soon after cyclic loading began, at 
approximately N=400 cycles. Very small delaminations were observed at the cracks, 
although, unlike specimens 2 and 3, these delaminations grew very slowly and stably with 
continued loading. In specimens 4 and 6, delaminations first grew toward the thick end of 
the laminate, and then grew from the crack toward the thin end. In specimen 5, the 
delamination grew only toward the thick end. The fatigue loading was continued for these 
laminates until a delamination had grown from the onset location to the grip at the fixed (thick) 
end. This condition was considered final failure. Figure 1 6 shows the final failure of 
specimens 4 and 5. After testing, the original photographs of specimens 4-6 were re- 
examined. Surface marcels were found at the locations where the initial cracks formed, as in 
specimens 2 and 3, although these marcels were of much smaller magnitude. 

Figure 17 compares the number of cycles to final delamination failure for specimens 2-6, 
with the results from ref. 2, for identical specimens, but without any ply waviness. The solid 
circular symbols in fig. 1 7 show the number of cycles to develop the initial crack through the 
wavy areas and the open circular symbols indicate the final failure. Even though the 
flexbeams in this study were tested at lower maximum strain levels than the flexbeams 
without waviness in ref. 2, failures occurred at significantly shorter lives. The fatigue lives of 
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specimens 2 and 3, with large amplitude marcels near the surface, were 1 to 2 orders of 
magnitude shorter than specimens 4-6, with smaller amplitude surface marcels. 

Analytical Results 

Ply waviness was characterized in this paper by the aspect ratio of the marcel, defined 
as the half-amplitude, a, divided by the width of the period, L (see the inset in fig. 18). In 
the FE model, the highest aspect ratios were found to be in the center ply-drop groups 
toward the thick-end of the flexbeam and in the isolated marcels near the surfaces. All of 
those marcels had an a/L of 0.055 to 0.085. The aspect ratios of the marcels in the center 
plies decreased gradually along the length of the flexbeam. 

Results of the FE model analysis were used to determine the peak interlaminar normal 
stresses, a n , in the model. The a n stresses reported are with reference to the element local 
(n-t) coordinate system, rather than the global coordinate system (fig. 10). Several studies 
have shown that delamination failure in composites with similar geometry is predominantly 
due to opening mode failure (normal to the ply direction) [2, 4,5, 6], with a much smaller 
contribution from the interlaminar shear mode. In this study, therefore, delamination was also 
assumed to be controlled by the interlaminar tension, and only the interlaminar normal 

stresses were considered. Peak a n stresses occurred in the model at marcel 1 , and another 

area of high a n stresses existed at marcel 3 (fig. 11). When the transverse load direction 
was reversed (V=4.45 kN, 1000 lbs), peak stresses occurred at marcel 2. In all cases, the 
peak stresses occurred at the most interior point of the marcel. 

In order to determine the effect of marcel location on interlaminar normal stresses, the 
model was inspected to find marcels of the same aspect ratio, but at different X- or Y- 
locations. Figure 18 shows the effect of varying the X-location of a marcel, while keeping 
the aspect ratio and distance from the surface constant. The figure compares calculated 
interlaminar normal stresses, at marcelled areas of the model with aspect ratios of either 
0.085 or 0.045, at approximately the same distance from the surface, but different X- 
locations. As the figure shows, for either aspect ratio, the stresses increase as the location 
of the marcel moves from the thick end toward the thin region. 

In fig. 19, the effect of the location of the marcel through-the-thickness is shown. For this 
case, stresses at different locations through-the-thickness are plotted for marcels with aspect 
ratios of 0.085 and 0.045, the same distances along the taper. The stresses are plotted as 
a function of the distance from the upper surface, y , divided by the section half-thickness, 
t(x). As fig. 19 shows, stresses decrease sharply for both aspect ratios, as the marcel 
location moves from the surface toward the midplane of the flexbeam. Exponential curves 
were fitted to the data in fig. 19. 

To determine the effect of the marcel aspect ratio, the aspect ratios of marcels 1-3 in the 
original FE model were modified and the analyses were repeated. The aspect ratio of 
marcel 1 was reduced from the original 0.085, to 0.055, 0.034, and 0.025. Marcel 2 was 
modified from the original aspect ratio of 0.05 to 0.028, 0.02, and then 0.014. The aspect 
ratio of marcel 3 was reduced from 0.085 to 0.08, 0.07, 0.06, and 0.03. Figure 20 shows a 

comparison of the a n stresses at locations marcel 1 , 2, and 3, as the aspect ratio was varied 
at each location. The peak interlaminar normal stresses increase exponentially as the aspect 
ratio increases. Again, the stresses are always highest for marcel 1 , the marcel furthest from 
the fixed end. The curves for marcels 2 and 3, which are near the thick end, and at similar X- 
locations, appear to cross at an aspect ratio of about 0.03, but at higher aspect ratios, 
marcel 3, nearest the thick end, is the lower bound of the stresses. 
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The test results from fig. 17 were used with the calculated curves from fig. 20 to correlate 
the interlaminar normal stresses with the fatigue lives of the flexbeam test specimens. For 
each of the tested specimens, the aspect ratio was measured at the marcel where the 
delamination failure initiated, the corresponding interlaminar stress from the FE analysis was 
determined from the appropriate curve in fig. 20. The curve for marcel 1 was used to 
determine the stresses for specimens 2 and 3, since for those specimens, delamination 
damage initiated there. For specimens 4-6, the curve for marcel 2 (at X=37.3 mm) in fig. 20 
was used, since the X-location was closest to the actual location of the marcels where failure 
began. For specimen 4, damage initiated at a small marcel at X=50.8 mm, and for 

specimen 5 and 6, at X=48.26 mm. The resulting a n values were plotted against the 
number of loading cycles to failure for the test specimens, in fig. 21 . As the interlaminar 
normal stresses decrease, the fatigue life increases. This figure also shows the data from 
the flexbeams without waviness from ref. 2. Finite element results from ref. 2 showed that 

the maximum value of a n under combined tension-bending loading was 12.5 MPa. Based 

on fig. 21 , a n = 13.0 MPa was chosen as the limit below which delamination failure will not 
occur at less than N=200,000 cycles. 

In fig. 22, this lower limit is plotted with the curves and data from fig. 20. The intersections 
of the curves show the maximum allowable a/L values for the corresponding X-locations in 
the flexbeam. Using the assumed lower limit, no marcel would be allowable at X=76.7. At 
the X=37.3 and X=27.4 locations, marcels with aspect ratios of no more than 0.015 and 
0.004, respectively, would be acceptable. 

Discussion 

Using this technique of combining test data with FE model results, a lower limit of 
interlaminar normal stress can be chosen for a desired fatigue life. That normal stress can 
then be used, as in fig. 22, to determine the maximum allowable aspect ratio at various 
locations in the flexbeam. 

The technique described above could be used as a first attempt at determining whether 
or not a given flexbeam should be used or discarded. Flowever, the method has several 
limitations: e. g., an unlimited number of wavy ply configurations are possible, but it is not 
practical to model each case to determine the effect on the interlaminar stresses. 
Furthermore, varying the mesh size in the FE model may change the calculated stresses 
significantly. A mesh refinement study is required to verify this. A more complete picture of 
the delamination process in these flexbeams may be provided by an energy-based 
analysis, such as that used in refs. 1 , 2, and 4, in which simulated delamination growth is 
modeled at various locations and the associated strain energy release rates are calculated. 
However, considering the discussion in fig. 17 and the considerably poorer performance of 
these flexbeams with marcels under the fatigue loading, the effort involved in such an 
analysis may not be justified. 

CONCLUDING REMARKS 

The effect of wavy plies on the durability of tapered laminated composite flexbeams 
was studied, for use in developing accept/reject criteria for these flexbeams. Tests were 
conducted on 25.4 mm wide coupon specimens, which were cut from a full-size flexbeam 
of S2/E7T1 glass/epoxy. All of the specimens had significant ply waviness through the 
center of the flexbeam and isolated areas with varying amounts of waviness at locations 
near the surfaces. The specimens were tested under combined axial tension and 
transverse bending loading. For all of the tested laminates, damage started as a crack 
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through a marcel near the surface. As the loading continued, delaminations grew from the 
crack in both directions along the length away from the marcel and multiple new 
delaminations formed at neighboring interfaces. Even though the flexbeams in this study 
were tested at lower maximum strain levels than identical specimens without ply waviness 
from ref. 2, failures occurred at significantly shorter lives. Fatigue lives for the tested 
laminates ranged from 10 3 to 10 4 5 cycles. In contrast, identical flexbeams without marcels 
had fatigue lives of 10 5 to 10 6 cycles at slightly higher bending loads. 

A parametric study was conducted to determine the effects of ply waviness along the 
taper and through-the-thickness. A (2-D) plane strain FE model was developed which 
modeled exactly the flexbeam lay-up and geometry of one of the test specimens, 
including the exact geometry of the wavy plies. The model was analyzed using the 
ABAQUS finite element program to determine the interlaminar normal stresses, a n . For 
marcels of the same aspect ratio, and at the same depth from the surface, a n stresses 
increased with distance along the taper, from thick end to thin end. For marcels of the same 
aspect ratio, and at the same location along the taper, the a n stresses decreased from the 
surface inward toward the flexbeam midplane. The model was modified several times by 
decreasing the amplitude of the waviness in three marcels near the surfaces and the 
analysis was repeated. It was found that the a n stresses increased exponentially as the 
aspect ratio was increased in these surface marcels. 

A technique was presented for determining the smallest acceptable marcel aspect ratio 
at various locations in the flexbeam. Test data were combined with calculated results for 
modeled marcels of different aspect ratios and at different locations, to determine a lower 
limit on the allowable interlaminar normal stresses. This limit was then related back to the 
calculated FE results to determine maximum allowable marcel aspect ratios to achieve a 
desired fatigue life. 
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Table 1. Material Properties 


Material 

E„, GPa 

E 2 2 , GPa 

G 12 , GPa 

V 12 

S2/E7T1 tape 

47.6 

12.6 

4.81 

0.28 

E-glass/E7Tl-2 fabric 

25.3 

24.1 

4.56 

0.153 

Steel 

201 

201 

77.3 

0.30 

• \ 2 ] 

Neat resin 

4.10 

4.10 

1.54 

0.33 


Table 2. Ply-Group Layups and Smeared Properties 


Ply Group 

Layup 

Exx, GPa 

Eyy, GPa 

G xy , GPa 

Vxy 

mid-plane 

[0 2 /45] s 

41.1 

13. 7 

6.32 

0.34 

Belt 1,2,4 

[0 4 ] 

47.6 

12.6 

4.81 

0.28 

Belt 3 

[±45/0 2 ] 

31.6 

15.6 

9.19 

0.33 

Fabric 

woven ±45 

25.3 

24.1 

4.56 

0.153 

Dropped 1-4 

[± 45 ] n 

14.8 

14.8 

13.3 

0.33 


Table 3. Strain Gage Locations 


Gage Number 

Distance from lower grip, mm (inches) 

1,6 

12.7 (0.5) 

2,7 

53.3 (2.1) 

3,8 

82.6 (3.25) 

4,9 

98.5 (3.88) 

5, 10 

160.0 (6.3) 
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Figure 4. Specimen 2 with wavy plies. 
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Figure 6. Axial tension and bending test stand and deformed flexbeam. 
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Figure 8. Maximum measured surface strain vs. applied transverse stroke. 
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Figure 11. Finite element mesh of tapered flexbeam with wavy plies and marcels near surfaces. 



strain gage # 



Figure 12. Maximum surface strains vs. flexbeam tip displacement. 





Figure 13. Surface strains in tapered flexbeam laminates under combined tension-bending loading. 
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Specimen 3, N =884 cycles 

Figure 15. Fatigue delamination damage in specimens 2 and 3 
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Figure 16. Fatigue delamination damage in specimens 4 and 5. 
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Figure 17. Fatigue delamination response of flexbeams with and without wavy plies. 












Figure 20. Interlaminar normal stresses at three locations with varying aspect ratio marcels. 





Figure 21. Calculated interlaminar normal stresses vs. cycles to delamination failure. 
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Figure 22. Marcel aspect ratio limits determined from test results and FE calculations. 
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